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Abstract 

Q . The Harper (or "almost Mathieu" ) equation plays an important role in stud- 

u. 

. ies of localization. Through a simple transformation, this equation can be 

converted into an iterative two dimensional skew-product mapping of the 
cylinder to itself. Localized states of the Harper system correspond to fractal 

> . 

, attractors with nonpositive maximal Lyapunov exponent in the dynamics of 

, the associated Harper map. We study this map and these strange nonchaotic 

attractors (SNAs) in detail in this paper. The spectral gaps of the Harper 
system have a unique labeling through a topological invariant of orbits of the 
! Harper map. This labeling associates an integer index with each gap, and the 

> ! 

. i-h , scaling properties of the width of the gaps as a function of potential strength, 

X. 

e depends on the index. SNAs occur in a large region in parameter space: 
these regions have a tongue-like shape and end on a Cantor set on the line 
e = 1 where the states are critically localized, and the spectrum is singular 
continuous. The SNAs of the Harper map are described in terms of their frac- 
tal properties, and the scaling behaviour of their power-spectra. These are 
created by unusual bifurcations and differ in many respects from SNAs that 
have hitherto been studied. The technique of studying a quantum eigenvalue 
problem in terms of the dynamics of an associated mapping can be applied to 



a number of related problems in 1 dimension. We discuss generalizations of 
the Harper potential as well as other quasiperiodic potentials in this context. 



I. INTRODUCTION 



The Harper equation, the discrete Schrodinger equation for a particle in a quasiperiodic 
potential, 



was originally introduced in the study of the motion of an electron in two dimensional lattice 
in a magnetic field Q. In the above equation, ip n is the wavefunction at lattice site n, uj 
is an irrational number, and 0o is a parameter. The nature of the eigenfunctions and the 
eigenspectrum has been of considerable interest since this equation arises in numerous and 
varied contexts 

Studies of localization phenomena have been important in condensed matter physics 
research for more than forty years ||, and here the Harper system has played a very im- 
portant role. In systems with a periodic potential, the Bloch theorem holds and leads to 
wavefunctions which typically are extended. In the contrasting case of a random potential, 
the states are localized. A quasiperiodic potential is intermediate between these limits, and 
it is of interest to determine whether states are localized or extended. Indeed, the problem 
of localization has important implications in a variety of fields, among the most significant 
being an explanation of transport properties in disordered solids f2p|| (for example the Mott 
theory of the metal-insulator transition f9j). In a 1 -dimensional tight-binding model with 
incommensurate modulation, at a critical strength of potential there exists a metal insulator 
transition. Such a quasiperiodic system is physically realizable, for instance, in quasicrystals 
which have incommensurate modulation and therefore lack translational order []1(J ) . 

As is evident, the Harper equation is invariant under the transformation ip n — > 
(— l) n i/) n , e — > — e, E — > —E, so that it suffices to consider the case E, e nonnegative. A 




(1) 
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number of studies P,|lCiHT7| have examined different properties of the eigenstates and spec- 
trum of the Harper equation. The most important feature of this system is that for e = 1, 



the equation is self-dual |L7|]: the equation for the Fourier transform of ip n is identical to 
Eq. flU). The quantum eigenstates are extended for e < 1, exponentially localized for e > 1 
and but at the critical value e = 1, the states are power-law localized. The eigenvalue 
spectrum is singular-continuous (namely, the spectrum forms a Cantor set) for all irrational 
values of u at this critical value of e M. For rational values of u, the Bloch theorem applies, 
and leads to a spectrum consisting of bands. 

In a study of the continuous version of Eq. ([[]), namely 

- d +2ecos(2irux)ip(x) = E^(x), (2) 
dx l 



Bondeson, Antonsen and Ott JI8] pointed out that by using Priifer's transformation, the 
above Schrodinger equation could be converted into an ordinary differential equation for the 
variable y = tan -1 (2gip'ifj) / ((tp') 2 — g 2 ip 2 ), g being an arbitrary constant, as 

% = - g {[g 2 + k 2 (x)} + [g 2 -k 2 (x)} C osy} (3) 

where k 2 (x) = E — 2e cos(27ro;x). This equation describes a forced dynamical system (the 
independent variable x is now the "time"), the forcing term deriving from the potential 
being quasiperiodic if u is irrational. Study of the orbits of Eq. @ gives information on the 
behaviour of wave-functions of Eq. (|2|). Indeed, the Lyapunov exponent A for an orbit of the 
dynamical system is related to the localization length 7 of the corresponding wavefunction 

0, 

A = -2/7. (4) 

Thus localised states(fmite 7 >0) correspond to attractors with negative Lyapunov expo- 
nents, while extended states (7 — > 00) correspond to orbits with zero Lyapunov exponent. 



As Bondeson et al. [1J| further argued, the quasiperiodic driving precluded the possibil- 
ity of a simple structure for these attractors, which were both fractal and had nonpositive 
Lyapunov exponents: these are strange nonchaotic attractors (SNAs) fT9f . 



Our concern in this paper is the dynamics of the Harper mapping which is a similar 
transformation of Eq. flU), using Ricatti variables, namely t/v-i/Vv - ► x n [^Oj^I| . Eq. flU) 
becomes the two-dimensional mapping 

x k+ i = -[x k - E + 2e cos 27r0 A .]~ i (5) 
4> k+ i = (f) k + OJ modi. (6) 

While this is exactly equivalent to the Harper equation, it is amenable to somewhat simpler 
analysis since boundary conditions that must be imposed for quantization in Eq. (|l|) become 
conditions on the dynamical state of orbits in the mapping, Eqs.(^-||). 

Viewed merely as an iterative mapping of the plane, the above system has several features 
of interest. The map is reversible, so there can be no chaotic dynamics, but for sufficiently 



large e, it can be argued [22] that there must be an attractor of the dynamics which has an 
infinite set of discontinuities. The attractors are therefore both strange (i.e. fractal) and 
nonchaotic (the Lyapunov exponents must be nonpositive): they are SNAs. The Harper 
map is another example of a quasiperiodically forced nonlinear dynamical system where it 
is known that SNAs are frequently encountered p3fl . By now a variety of quasiperiodically 
driven systems with SNA dynamics are known; these include examples of driven pendu- 
lumsce, forced nonlinear maps, the driven van der Pol oscillator, neural network systems, 



as well as some mappings where there is no explicit quasiperiodic driving [23]. Owing to 



the fractal structure of the attractors, the dynamics is strictly aperiodic, but the correlation 



properties p4J of SNAs show them to be much less aperiodic than chaotic attractors, the 
correlation function always lying between and 1. 

Strange nonchaotic dynamics is of considerable interest since it is intermediate in several 
ways between regular and chaotic motion. SNAs are typically found in the region of a 



transition to chaos, and several routes or scenarios for their formation are now known |23 
Their existence can be confirmed by explicit calculation of Lyapunov exponents and fractal 
dimensions, or by computing their spectral and correlation properties. In the Harper system 
it is known that SNAs can be created by the homoclinic collision of an invariant curve with 



itself | 25| ) and that this route to SNA is accompanied by an unusual symmetry breaking. 
There are other routes through which SNAs can also be formed, and some of these are 
described here as well. 

The Harper system is therefore of interest from both the nonlinear dynamics perspective, 



as well as in the context of localization. Indeed, as initial studies have shown p2| , |25| , |26 



both these aspects of the problem are related, and the correspondence between SNAs and 
localized states can be exploited so as to provide new insights in the study of forced nonlinear 
systems, as well as for localization studies. It has been noted for instance, that quasiperi- 
odic potentials are frequently associated with critical or power-law localized states, just as 
quasiperiodic driving frequently transforms chaotic attractors to SNAs. This can be made 
more concrete by studying the iterative mappings analogous to Eqs. that are associated 
with such potentials. In many instances, it also turns out that such systems have attractors 
which are SNAs [|7|. 

This paper is concerned with such correspondences. One focus is on what can be learned 
about the spectrum and eigenstates of the Harper equation through a study of the dynamics 
of the Harper map. At the duality point, e = 1, the eigenspectrum is given exactly by the 
zeros of the Lyapunov exponent, \(E). Furthermore, topological properties of the orbits of 
the mapping provide a unique labeling of the gaps in the spectrum. The gap indices also 



determine the scaling behaviour of the g function of the coupling constant p6[ . These 

results are described in Section III of this paper. The behaviour of the gaps as a function 
of coupling constant is important for a number of problems, and we have shown [2B] how 
the gaps can be uniquely labeled by mapping them on to a Cayley tree. The gap labels are 
winding numbers, which are topological invariants for orbits of the iterative mapping. These 
determine the scaling behaviour of the gaps at, above, and by duality, below, the critical 
coupling, e = e c = 1. We discuss the general scheme for labeling the gaps for particular 
choices of lu such as the golden mean or silver mean ratios, as well as for other algebraic and 
transcendental irrational ratios. 
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At the same time, the SNAs of the Harper map are themselves of interest. Regions 
(in parameter space) where they occur have an interesting hierarchical geometry, being 
organized in fractal "tongue-like" zones [p6fl . In Section II of this paper, we first describe 
properties of the iterative mapping for the Harper potential. The phase diagram for this 
system, as a function of E and e for fixed irrational uj is discussed. The different bifurcation 
routes for the formation of SNAs, and the measures have been introduced in the past few 
years in order to distinguish SNAs from other related attractors are elaborated in Section 
IV. 

The case of quasiperiodic potentials obtained by generalizing the Harper potential, where 
critical states are known to occur, are discussed in Section V. We show that the technique of 
determining the spectrum from the behaviour of the Lyapunov exponent versus parameter 
curve can be extended to these potentials as well. Other cases of systems with critical states 
are known, as for example the quasiperiodic potentials that arise from the Fibonacci, Thue- 
Morse, Rudin-Shapiro and period-doubling symbolic sequences. The corresponding iterative 
mappings also support SNAs f28fl. The paper concludes with a summary in Section VI. 



II. THE HARPER MAP 

We study the mapping Eqs. (|5|-§|) and other quasiperiodic generalizations of it which 
arise from quasiperiodic Schrodinger problems. While most of the essential features of the 
system do not depend on the precise value of uj so long as it is irrational, some of the details 
(such as the arrangements of the gaps, see Section [Hj below) do. In this work uj is taken 
to be a quadratic irrational number such as the golden mean ratio — l)/2, which is 
the value for which most previous studies have been carried out, or the silver mean ratio, 
\/2- 1. 

The Harper map is a skew-product dynamical system on the cylinder S 1 ® R. Since the 
4> dynamics is quasiperiodic, there are no periodic orbits. There is also no chaotic dynamics. 
In general there are three qualitatively different types of motion that can occur in this 
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system. There can be quasiperiodic tori (two-frequency or three-frequency) with negative 
Lyapunov exponent, neutrally stable curves with zero Lyapunov exponent, when the phase 
space is conserved under the flow, and strange nonchaotic dynamics. On SNAs the motion 
is both aperiodic and, as is characteristic of such attractors, intermittent 



The <fi- freedom is neutral with the corresponding Lyapunov exponent being trivially 
zero, while for the x- freedom, the stretch exponent or one-step Lyapunov exponent BO] is 



yi = ln|x- +1 | = 21n|x i+ i|. (7) 
The iV-step finite-time Lyapunov exponent is 

^ N-l 

*N{Xi) = JaYl VjM- ( 8 ) 
3=0 

These local exponents are indexed by initial condition x iy but the asymptotic Lyapunov 
exponent, 

A = lim X N (xi), (9) 

iV^oo 

is independent of (almost every) initial condition. In the rest of the paper, we deal only 
with this nontrivial Lyapunov exponent, for which the following relation holds fTWj 



X(E,e) = — 21n(e) + X(l/e,E/e). (10) 

relating the exponents above and below the critical coupling. 

The eigenvalue E appears simply as another parameter in the Harper map. Quantization 
conditions that apply for the Harper equation, namely that the wavefunctions should be 
normalizable, now become conditions on the dynamical state of the map. In particular, 
these can be viewed as conditions on the Lyapunov exponents: for extended states the 
Lyapunov exponent is zero, while for localized states, the Lyapunov exponent is negative. 

At the duality point e c , eigenstates of the Harper equation are critically localized. Hence, 
the Lyapunov exponent of the Harper map is identically zero whenever E is an eigenvalue. 
Furthermore, since the spectrum is known to be singular-continuous, the A versus E curve 



meets the line A = on a Cantor set of points (Figs. la,b). Indeed, the Hofstadter butterfly 
@] which is the spectrum of the Harper equation as a function of uj can be determined as 
the set {E : X(E,uj,e c ) = 0}. 

The variation of the Lyapunov exponent with e is simple if E is an eigenvalue. Below 
e = 1, the exponent is zero and the state is extended, while above it, A ~ — lne and the 
state is localized; see Fig. 2 for the case of E = 0. Of particular note is the discontinuity 
(in slope) at e c : the localization transition can therefore be viewed as a bifurcation. As 
shown previously, this bifurcation is accompanied by a quasiperiodic symmetry-breaking 
and the homoclinic collision of neutrally stable invariant curves to form a SNA ||25|| . For 
other eigenvalues, the location of the eigenvalue changes with e. At e c there is always a 
symmetry-breaking, but there need not be the homoclinic collision, so SNAs can be formed 
here by other routes [[H],^ as is discussed in Section [TV] below. 

The overall dynamical behaviour in the Harper map can be summarized in a phase 
diagram, shown in Fig. 3a for the case of uj = y2 — 1, namely the silver mean ratio. Owing 
to the E — ► —E symmetry, it suffices to consider (E > 0, e > 0). We compute the Lyapunov 
exponent and determine the dynamical state in the different regions demarcated in the 
figures. Below e c , the dynamics is either on invariant curves (I) or on tori (Q), while above 
e c , there are regions of SNA (S) and tori (Q). The SNA regions exist only at and above e c . 
As was first argued by Ketoja and Satija P2"l , in the strong-coupling limit of large e, the 
dynamics is guaranteed to be on SNAs; this argument can be extended down to e c . Above 
e c , therefore, the dynamics is largely on SNAs. The gap regions extend above e c . Therefore, 
by continuity, the SNA regions emanate from the eigenvalues at e c . These have a tongue-like 
shape, and since each eigenvalue is associated with one such region, these are hierarchically 
organized in the same manner as are the energy eigenvalues: SNAs occur in fractal tongues. 
Shown in Fig. 3b is a (schematic) detail in the neighborhood of one (and therefore of every) 
tongue: the region of 3-frequency quasiperiodic orbits or tori for e < 1 are mirrored by SNA 
regions for e > 1. 
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The major differences between the cases of different u arise from the manner in which 
the regions of tori (Q) are arranged. These correspond to gaps in the eigenspectrum of the 
Harper operator, which we now discuss. 



III. GAPS OF THE HARPER SPECTRUM 

The eigenvalue spectrum of the Harper equation is singular-continuous for irrational 
u at e c , namely it forms a Cantor set. The nature of the gaps in the spectrum has been of 
great interest in a number of different contexts in the past several years. In particular, how 
the gaps behave as a function of the potential strength, namely e, has been a long-standing 
question p0| , |33 |. 



The eigenvalues can be simply deduced from the A versus E curve at e c , as shown in 
Fig. la for the case of u the golden mean ratio, or in Fig. lb for w the silver mean ratio. 
The tori regions of the Harper map correspond to gaps in the eigenvalue spectrum of the 
Harper equation. Orbits for such values of E and e are 1-d curves that wind across the 
cylinder an integral number of times. Within each gap, this integer index is constant and 
measures the number of times the orbit crosses the boundaries x — > ±00. This essentially 
counts the number of changes of sign of the wave-function (per unit length), and is related 
to a winding number for the orbit. These indices are indicated in Figs, la and lb for the 
largest visible gaps. 

For rational frequencies, the spectrum of the Harper equation consists of a finite number 
of bands. If to = p/q,p, q integers, then there are q or q — 1 bands depending on whether q is 
odd or even, and therefore in the positive energy spectrum the number of gaps is (g— 1)/2 or 
(q/2 — 1). For rational uo all orbits are periodic with period q. For energies in the gaps, the 
orbits form a 1-parameter family, indexed by initial phase, from which the winding number 
can be easily extracted. 

As is well-known pi}] , an irrational number has a unique continued fraction representa- 
tion. Thus an irrational frequency 
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Gtj's integers, can be successively approximated by a unique series of rational approximants. 
By truncating the continued fraction representation up to the term a,j, one obtains the 
rational Pj/qf, the sequence Pj/qj,j = 1,2,... converges, Um^ooPk/qk — ► w. For a given 
approximant Pj/qj, the gap structure of the spectrum can be described in terms of the 
winding numbers. These run from 1 to qj/2 — 1 or (qj — l)/2, as discussed in the preceding 
paragraph. The new gaps that appear for the next approximant, pj+i/qj+i, do not alter the 
inter se ordering of the gaps already present for lower order rational approximants. Thus 
the gap structure for irrational u can be deduced by sequentially examining the structure 
of the gaps for the successive rational approximations to uj. Each gap in the spectrum, or 
each region of tori in the map, can therefore be uniquely labeled by an integer index. 

For algebraic irrationals which have periodic continued fraction representations, there 
are simplifications that permit a more compact description of the ordering scheme. We now 
describe this for the case for the golden and silver mean ratios. 

The golden mean is the limiting ratio of successive Fibonacci numbers, Fk/Fk+i, where 
Fk is given by the relation 



fc+i 



F k + F k -!, k = 1, 2, ... , with F = 0, F 1 = 1, (12) 



the first few being 0,1,1,2,3,5,8,. . . , while the silver mean is the limiting ratio of the successive 
integers of the family 0,1,2,5,12,29,70,. . . , namely 

S k+1 = 2S k + S k -i, k = l,2,..., with Si, = 0, S x = 1 (13) 

with Hindoo Skj Sk+i — >• \/2 - 1. 

The ordering of the gaps can be best described by first constructing a graph for the 
particular irrational ratio in the following manner. We discuss this in detail for the silver 



mean ratio below, and indicate the analogous result for the golden mean ratio [26 . 
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Consider a Cayley tree, arranged as shown in Fig. 4a, with each node (except the origin, 
labeled 0) having two successors. Nodes at the same horizontal level are at the same gener- 
ation. The rightmost node at each generation is labeled by successive <S&, while the leftmost 
are half the successive even <S&'s. The other integers are identified with the remaining nodes 
as follows: for given m, the mother node i m is the smallest available such that 

• m + i m is one of the integers Sj.- If im already has two daughters, then the mother 
node is chosen such that 

• m + i m is one of the integers Sk+i = Sk+ 2 — <Sfc+i, namely 1,3,7,17,41,. . . . 

For u the golden mean, i m is chosen such that m + i m is the smallest possible Fibonacci 
integer [25]. Complicated (but deducible |25fl ) rules determine which branch (left or right) 



the daughters are placed, but the above procedure results in an unique arrangement of all 
integers on the Cayley tree: see Figs. 4a and 4b. 

The actual ordering of the gaps derives from the following additional consideration. Every 
pair of integers, i\ and with i 2 > h, has two possibilities as to how they are relatively 
placed on this graph. Either 

1. i\ is an ancestor of i 2 , i.e. there is a directed path connecting i 2 to i\. If this path is 
to the left at node i\, then i 2 -< i\. (If to the right, then i\ -< 

or 

2. z is the most recent common ancestor of i\ and i 2 - If the path from z to i\ is on the 
left at io, then i\ -< i 2 . (Similarly, if it is to the right, then i 2 -< 

This gives a unique ordering of the integers with the relation -< being transitive (if i -< j 
and j -< m then i -< m) . Thus, for the silver mean we get 

...-< 35 -<...-< 11 -<... -<l-<...-<4-<...-<2-<...-<7-<...-< 12 -<...-<0. 

while for the golden mean, the ordering is 

...^4^...^9^...^1^...^7^...^2^...^11^...^21^...^0. 
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The gaps appear in precisely this order: if k -< £, then the gap with index k precedes 
the gap with index t in the positive energy spectrum of the Harper system: see Figs. la,b 
where gap labels for the largest visible gaps are indicated. The labels have been determined 
by examination of orbits of the map as described. The attractor shown in Fig. 5 is an orbit 
in the gap labeled 8 (which lies between gaps 4 and 3) in the silver-mean Harper system. 

As can be seen in Figs, la and lb, each gap has two attributes: the width and the depth. 
The width is the spacing between two energy eigenvalues of the Harper equation, but the 
depth has no immediate significance, being just the minimum value that the Lyapunov 
exponent takes between two zeros. These are denoted w m and d m respectively, and both are 
functions of the coupling strength, e. 

In general, both widths and depths decrease with order (i.e. with increasing gap index), 
and with increasing e. However, and particularly in the case of the widths, the dependence 
on gap label is nonmonotonic, in a manner described below. The depths scale in a simple 
fashion, namely 

d » ~ ]vvv> ( 14 ) 
as shown in Figs. 6a,b. On the other hand, the widths must be considered in families labeled 
along a consistent branch of the Cayley tree. For example, for uj = v2 — 1 (see Fig. 4a), 
the gaps 1,2,5,12,29,. . . , the rightmost nodes at each generation form one family, while gaps 
1,6,35,204,. . . , the leftmost nodes, form another family, for each of which the widths scale, 
at 6^ , as 

w N ~ ^, (15) 

the exponent 9 being particular to a given family of gaps. The two families of gaps enumer- 
ated above are the fastest and slowest decreasing, respectively, with exponents 9i ~ 2.05 and 
9 r ~ 2.31. For all other families, the exponents lie within the limits [9i,9 r ]. Similar results 
for the golden mean case are shown in Fig. 5b. 

Wiegmann et al. |§ have studied the tight-binding Harper Hamiltonian for which they 
obtain the gap distribution, namely the number of gaps lying between D and D + dD. They 
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find that this quantity, p{D) ~ D 7 where 7 is found to be approximately 3/2; similar 



exponents have been reported by Geisel et al. [15[ as well. If the gaps are rank-ordered, 
namely from largest downwards, disregarding the gap indices, then we observe that they 
decrease as 

w r ~ 1. (16) 
// ps 2, which is consistent with the power-laws obtained earlier P,|T5|]. 



IV. SNAS IN THE HARPER SYSTEM 

We now turn to the description of the dynamical attractors in the Harper map. Strange 
nonchaotic dynamics in the Harper system is known to be created through a number of 
different routes or scenarios. Unlike other examples of quasiperiodically forced systems 
wherein SNAs are formed by mechanisms that have parallels with routes to chaos in similar 
unforced systems, here the routes to SNA can be quite distinctive. 



A. Homoclinic collision and symmetry breaking 

As has been described earlier, at band center, namely for E = 0, orbits of the Harper 
map below e c lie on invariant curves that consist of two branches (these are denoted C for 
central and N for noncentral in Fig. 7a). These branches originate in the fact that the 
mapping for e = has a period-2 fixed point and no period-1 fixed point. Different initial 
points on the (x, <fi) plane generate different invariant curves, the two branches of which 
are widely separated at e=0, but as e — * 1, these branches begin to approach each other. 
For a particular curve, the branches collide at a dense set of points, at e c , and indeed, 
all the different curves merge at e c to give a single attractor which has a dense set of 
singularities (Fig. 7b). At the merging point, the Lyapunov exponent remains zero. The 
distance between the two branches of a given invariant curve decreases as a power in e c -e, as 
does the distance between any pair of invariant curves (Fig. 7c). As previously noted, the 
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homoclinic collision route to SNA is accompanied by a symmetry-breaking. Below e c when 
the Lyapunov exponent is zero, there is a quasiperiodic symmetry in the stretch exponents 
(Fig. 7d), and this symmetry is lost above e c , once the attractor has negative Lyapunov 
exponent (Fig. 7e). 

While each eigenvalue at e c is associated with a SNA, not all SNAs arise from the above 
homoclinic collision route. Indeed, the more general route to SNA in this system appears 
to be related to both the blowout bifurcation as well as the so-called fractalization scenario 
for the formation of SNAs, and this is discussed in the following subsection. 



B. Other Bifurcations and Fractalization 

Given the complexity of the phase diagram for the Harper system with the hierarchical 
organization of SNA (S) and torus (Q,G) regions, there are a number of bifurcations as 
parameters are varied. For fixed e below e c , there are bifurcations from 2-frequency to 
3-frequency tori as E is varied, the signature of this transition being the change of the 
Lyapunov exponent from a negative value to zero. There is therefore a discontinuity in the 
slope of the Lyapunov exponent versus E curve. The 3-frequency tori (Q) densely cover 
the phase space, while the 2-frequency tori (G) are attractors of the dynamics (including 
the case when the Lyapunov exponent is exactly zero). This transition is also, therefore, 
accompanied by a breaking of the symmetry in the stretch exponents. On the other hand, 
by fixing E and varying e, there are other bifurcations depending on the value of E. There 
can be transitions from 3-frequency tori to 2-frequency tori below e c , and from 2-frequency 
tori to SNA above e c . 

Shown in Fig. 8a is the variation of the Lyapunov exponent with e for E — 1; the 
bifurcations being signaled by a discontinuity in slope of the curve, with the Lyapunov 
exponent being strictly nonpositive throughout. SNAs are first formed at this energy at 
e~ 2.5200, below which there is an attracting 2-frequency torus. This transition has some 



of the dynamical characteristics of the blowout bifurcation |35| both in the nature of the 
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attractors before (Fig. 8c) and after (Fig. 8d) the transition as well as in the manner in 
which the Lyapunov exponent varies. The fluctuations in the local Lyapunov exponents, as 
measured in the variance of the distribution, are large on the SNAs and zero on the tori. 
This is shown in Fig. 8b. Note that the Cantor set structure of the gaps (whose widths 
decay only as a power in e) is reflected in the pattern of bifurcations along any typical line 
in the (E,e) plane. 

The G— > S transition occurs along essentially all the boundaries of the SNA regions above 
e c . When approached from below, the attractor of the dynamics can be seen to gradually 
develop from a smooth curve to one with wrinkles on all scales: this is reminiscent of the 
fractalization process [fj^l . The blowout nature of the transit ion-insofar as the attractor 
abruptly changes its volume with intermittent "bursts" -can be seen more clearly along the 
line E = 3. In the torus region G, the dynamical map here has a single invariant curve 
(Fig. 9a), which develops kinks and wrinkles as e is varied. The transition to SNA is at 
e ~ 1.270(37567), when the attractor gets wrinkled enough and actually hits the boundary 
x — > oo, Fig. 9b. This bifurcation can be clearly detected by measuring, say, the distance 
A = min„ |1 — tanhx n |, namely the closest point of the attractor from the boundary as a 
function of e. As can be seen in Fig. 9c, this quantity abruptly decreases at the transition. 



C. Characterization of Harper SNAs: Correlations 

Principal among the measures that can be used to distinguish SNAs from morphologically 
similar chaotic attractors or from (regular) attractors which are similar dynamically (namely 
they have similar values of the Lyapunov exponent) are quantities that examine correlations 



on the attractors 24,36 



The autocorrelation function, which provides a quantitative measure of the extent to 
which dynamical properties are correlated, is defined for a dynamical variable x as 

C(T) = «> - <«.)' (17) 
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where i = 1, 2, ... is a discrete time index, r = 0, 1, ... is the time shift, and () denotes a 
time-average. If x is periodic then C{r) varies between -1 and 1, regaining its initial value 
C(0) in a periodic fashion. If x is chaotic, then C(r) decays from its initial value of 1 to 0, 
around which latter value it fluctuates. On SNAs, the variables have correlation functions 
that do not recur exactly: C(r) is also quasiperiodic and is nearly recurrent, oscillating 
between -1 and 1 (Fig. 10a). 

The Fourier transform of the correlation function also shows these differences as charac- 
teristic spectra. The Fourier transform of a discrete sequence {xk} is 



1 N 

T(Q,N) = -5> fc exp(z27rK]), (18) 
fc=i 



the power spectrum being 



P(w)= lim(|T(^,iV)| 2 ). (19) 

N^co 



A chaotic signal typically has a continuous power spectrum, but for SNAs, the Fourier 
spectrum reflects the fact that the dynamics is neither chaotic nor regular, and the power 
spectrum consists of several peaks at frequencies corresponding to the quasiperiodicity in 
the dynamics; see Fig. 10b. 

D. Phase Sensitivity Properties 



Pikovsky and Feudel |24| introduced the phase-sensitivity exponent in order to charac- 
terize the strangeness of a attractor. This quantity measures the sensitivity of the dynamics 
with respect to change of phase of the external force. Given a map 

■En+l f \%n ; @n) 

9 n+1 =9 n + co mod 1, (20) 

the maximal value of the derivative dx/89 can provide a suitable tool to distinguish between 
strange and nonstrange geometry. Since the above derivative changes along the trajectory 
(xq, 6*o), (xi, 9i), (x%, 62), ■ ■ ., one can get a relation 
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o dx n+1 

= f(x n ,e n ) + f( Xn ,e n )^- (21) 

where Sn is the derivative with respect to external phase (we follow the notation of Pikovsky 
and Feudel |24|). The number of maxima in Sn up to time N is 

1 N (x,0)=max\S n \. (22) 

0<n<iV 

The value of 7^ grows with N, a consequence of the fact that the attractor is not smooth; 
see Fig 11a. Also, as Pikovsky and Feudel |24| noted, the growth rate of the partial sum 
with time represent the strangeness of attractor, so it becomes necessary to calculate the 
minimum value of 7at(^, 9) with respect to randomly choose initial points, 

T N = mm[ max :\S n \ ]. (23) 

x,0 0<n<N 

From Fig. lib we see that Tn ~ N^, where /i is the so-called phase sensitivity exponent. 
If the attractor is smooth (i.e. nonfractal) then the maximum derivative with respect to 
external phase is bounded, and this value saturates as a function of iteration. However, in the 
case of SNAs, the maximum derivative with respect to external phase increases indefinitely 
with iteration. 



V. MAPS OF HARPER TYPE 

In this section we discuss the dynamics of mappings similar to Eqs. (HH). The Harper 
map can be generalized in two essentially different ways. By generalizing the potential in the 
Harper equation but retaining the feature of quasiperiodicity, one essentially obtains other 
Schrodinger equations. Power-law localization is known to occur in a variety of quasiperiodic 
systems, well studied examples being the potentials J3?[ 

V n — ei cos 2irnuj + e 3 cos 67rnuj (24) 

and 
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V n = — tanh(ucos27mc<A (25) 

tanh \i 

In both these systems, it is known that there are mobility edges. In the resulting maps, 
which are very similar to the Harper, the cosine "driving" term is replaced by a more 
general function which is also quasiperiodic since u is irrational. 

Since both these above cases can be viewed as modifications of the Harper system, the 
phase diagrams obtained by variation of E and the parameters (ei, €■$ or e, fi) are quite similar 
to Fig. 3a,b. Some instances have been studied previously. For Eq. ( |24"D with e 3 = 1/4, the 
regions of 3- and 2-frequency tori are arranged in a manner very similar to the case with 
€3 = 0, namely the Harper map (Fig. 3a), except that the transition from extended to 
localized states does not occur on the line ei = 1, but instead on some complicated curve 
in the (E, e\) plane; see Fig. 12a. Fig. 12b shows the corresponding phase diagram for the 
potential in Eq. ([25|). 

A number of quasiperiodic potentials generated from abstract sequences support critical 
states [p0|| . Among the most extensively studied such sequences is the Fibonacci chain, the 
potential for which is defined by 

V n = a < {nuj} < uj 

= -a u < {nu} < 1. (26) 

For any value of a, all states of this system are critical [P^-HD[. The corresponding 



Schrodinger equation can be transformed to the mapping 

x k+1 = -[xk-E + Vk]- 1 (27) 

with Vk given by Eq. (^). Analogous maps can be obtained for the Thue-Morse, Rudin- 
Shapiro and period-doubling quasiperiodic sequences, and these have been shown to support 
SNAs [P7[ . The nature of the eigenspectrum in such cases will be discussed elsewhere . 
Another possible generalization of the Harper map gives 

Xk+i = -\fM + 2ecos27r0 fc ]- 1 (28) 
4> k+ i = 4> k + uj mod 1. (29) 



where f(x) is now arbitrary while the quasiperiodic driving term is retained. Because of 
the somewhat unusual nature of the mapping, not much can be said in general — there can 
now be chaotic motion as well if / is not 1—1, for example. However, analysis of the strong- 
coupling limit still suggests that SNAs should arise in such systems. Indeed, for special cases 



that have been examined so far p5| , |3T| many of the scenarios to SNA that have been seen in 
the Harper mapping also exist in these more general cases. An important limitation is that 
except for the case when / is linear, the above mapping cannot arise from a Schrodinger 
equation, and thus the parallel insights that derive from study of the quantum-mechanical 
system are not available for nonlinear f(x). 

VI. SUMMARY AND DISCUSSION 

In the present work we have studied the Harper equation from a dynamical systems point 
of view. The correspondence between localized states and strange nonchaotic attractors 
jj|22] makes it possible to obtain a number of properties of the eigenspectrum, as well as 



the eigenstates of the Harper equation through a study of the equivalent iterative Harper 
map. 

For the spectral gaps, we have described a labeling scheme which explicitly connects 
the number-theoretic properties of u with the gap ordering. This is most conveniently 
depicted on a Cayley tree, which we describe in detail for the case of the so-called silver- 
mean quadratic irrational cu = \[2 — 1. The general principles governing this organization 
are not difficult to ascertain, and computationally, these are most easily deduced through 
examination of the orbits of the Harper map. The gap indices that we obtain are related 



to the winding numbers introduced earlier by Johnson and Moser fril l, an d these (integer) 
indices determine the scaling behaviour of the g function of the coupling. We find 

that the width of each gap decreases only as a power of the coupling, and thus the gaps 
never close: this confirms the conjecture of Thouless et al. ||33|| . The organization of the 
gaps on the Cayley tree also makes it possible to identify other scaling properties. For fixed 
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e, the widths of gaps along a branch of the tree also show a power-law scaling in the gap 
indices. 

The transition from extended to localized states is accompanied by a breaking of the 



symmetry of stretch exponents \ 25[] . This symmetry-breaking appears to accompany a num- 
ber of the bifurcations in this system. Recall that the stretch exponents are the logarithms 
of the ratio of wavefunctions at adjacent sites. Thus the symmetry that gives a zero value of 
the Lyapunov exponent in the extended state is the following: the ratio of the wavefunctions 
at a triple of sites, i — + (a:b:c, say) is such that at the triple of sites j — 1, j, j + 1, the 
ratio of wavefunctions is (approximately, because of the quasiperiodicity) reversed (namely 
c:b:a), for each j such that j — i is a Fibonacci number. On the other hand, the wave- 
functions of the localized phase do not have this symmetry and are instead characterized 
through their multifractal oscillations, which have also been extensively studied |[42|| . 

There have been several realizations of the Harper system in experiments. Experiments 
of wave-chaos in microwave scattering [43[] in a geometry that is designed to mimic the 



Harper equation have shown the Hofstadter butterfly. In a very recent experimental study 



of the Hall conductance of a suitably designed heterostructure |44| , the Cantor-set nature of 
the energy spectrum has been explicitly detected. The scaling results for the gaps reported 
in this paper and elucidation of the gap organization are relevant to such studies. 

The methods of analysis employed here have wide applicability. A number of quasiperi- 
odic systems in 1-dimension are known to have singular continuous spectra. The correspon- 
dences studied here for the Harper system appear to be valid in the general cases as well, 
and thus offer an alternate and powerful means of analyzing such problems. 
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Figure Captions 

Fig. 1. (a) Lyapunov exponent versus energy at e c for uo = — 1)/2. Gap labels are indicated 
for the largest visible gaps. At every bifurcation, when A = 0, the dynamics is on a 
SNA. (b) Lyapunov exponent versus energy at e c for uo = \[2 — 1. Gap labels are 
indicated, as in Fig. la, for the largest visible gaps. 

Fig. 2. Lyapunov exponent versus e for the eigenvalue E — 0. The extended to localized 
transition occurs at e—e c —l. 

Fig. 3. (a) Phase diagram for the Harper map for uo = \[2 — 1. Below e c , namely the vertical 
line the dynamics is on three-frequency quasiperiodic (Q) orbits or extended states, 
1-d attractors which are two-frequency quasiperiodic orbits or gaps (G), and above 
e c , the dynamics is on SNAs (S), and 1-d attractors (G). Only the largest gaps are 
visible at this scale. 

(b) Detail in the neighborhood of the central gaps. Every tongue of SNAs necessarily 
has a similar structure. 

Fig. 4. (a) Ordering of the gaps for uo = (\/5 — l)/2, the golden mean ratio. Only part of the 
Cayley tree described in the text is shown for clarity. Each node has two daughters 
except for 0, which has only one. 

(b) The corresponding Cayley tree for the case of uo the silver mean, \[2 — I. 

Fig. 5. The attractor for a value of E corresponding to the gap N = 8. Note that the orbit 
has 8 branches that traverse the range — oo < x < oo. 

Fig. 6. Scaling of the gap widths, w N (•), and depths d N (o) as a function of gap index, N, 
at e = e c for the cases of (a) uo = (y/b — l)/2 and (b) uo = \[2 — I. For clarity, the 
depths have been multiplied by a factor of 10. The dashed line fitting the depths has 
slope -1. The dotted lines show the scaling of the two families of gaps; see the text for 
details. 



Fig. 7. (a) Invariant curve for the Harper map corresponding to the band center extended 
state, the central and noncentral branches being indicated by C and N respectively, 
(b) The fractal attractor corresponding to the critically localized state, (c) Scaling of 
the minimum vertical distance d between a family of different invariant curves in the 
Harper map as a function of the parameter (1 — e) for the eigenvalue E = 0. Different 
pairs of curves approach each other at different rates, but at e c , all curves collide 
to form the critical SNA. (d) First return map for the stretch exponent showing the 
quasiperiodic symmetry below e c , and (e) the return map above e c , when this symmetry 
is broken. 

Fig. 8. The depressed blowout bifurcation along the line E = 1. (a) Variation of the Lyapunov 
exponent with e for u the golden mean, (b) Fluctuations of the Lyapunov exponent. 
The attractor (c) before (a two-frequency torus) and (d) after (a SNA) the bifurcation. 

Fig. 9. The bifurcations along the line E — 3, showing the attractor (a) before (e = 
1.270375 . . .) and (b) after the G —>■ S transition, (c) The distance from the attractor 
to the boundary function of e. 

Fig. 10. For the E = critical SNA for uj the golden mean (a) the correlation function, and 
(b) the power spectrum. 

Fig. 11. (a) Maxima of the partial sum 7^ (cf. Eq. (JZ2])) on the critical SNA at e c , E=0 in the 
Harper map, as a function of iteration length, N. (b) T^, the minimum value of 7^. 

Fig. 12. Phase diagram for (a) the generalized Harper model, Eq. ([24]) with e\ = and €3. The 
tongue-like regions where SNAs can be found continue to be hierarchically arranged 
but the transition from localized to extended states occurs along a curve in the plane 
rather than on the line e=l. (b) As in Fig. 12(a) for the potential of Eq. fl25|) . 
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